Theory of Coupled Phase Transitions: Phase Separation and Abnormal Variation of 

Order Parameter 
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A simplified Ginzburg-Landau theory is presented to study generally a coupling of a first-order phase 
transition (FOPT) to a second-order phase transition (SOPT). We show analytically that, due to the 
coupling between the two phase transitions, the SOPT may exhibit a FOPT-like phase separation 
in which an ordered phase is separated from a disordered one. This phase separation results in a 
distinct behavior in the variation of the order parameter of the SOPT, namely, it is primarily the 
proportion of the ordered phase that contributes to the total order of the whole system during the 
FOPT. This growth mode may turn a mean- field critical exponent from 1/2 to 1 or even bigger. 
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Phase transitions have been one of the fascinating 
fields in statistical and condensed matter physics, in 
which a recent theme is the behavior of transition metal 
oxides akin to the cuprate superconductor, a topic which 
fills most part of recent journals in association with such 
rich phenomena as charge and spin orders, stripe phases 
and phase separations. In complex systems, various 
phase transitions have been detected experimentally to 
coexist such as difi^usion and shear ||l[ , superconductivity 
and ferromagnetism |^ , and antiferromagnetism and fer- 
roelectrics 0. Meanwhile, electronic or magnetic tran- 
sitions may compete with structural changes [Q. For 
systems containing polymers like monolayers and liquid 
crystals, positional and orientational orders frequently 
come into play As more sophisticated experimental 
techniques become available, more complex phenomena 
will certainly be revealed with more details, a theoreti- 
cal understanding of such complex phenomena within the 
theory of phase transitions is therefore highly demanding. 

In a recent Letter we put forward a Ginzburg- 
Landau theory based explicitly on the concept of cou- 
pled phase transitions. Taking the bainitic transforma- 
tion as an example, we showed numerically that the cou- 
pling of two FOPTs representing respectively the diffu- 
sional and displacive transitions in a carbon-ion system 
leads to such possible phenomena as diffusional-phase- 
transition-induced nucleation of the displacive transition 
and incomplete reaction. Here, employing a simplified 
free-energy functional, we are able to show analytically 
that for a coupling of a FOPT with a SOPT, the latter 
transition may display a phase separation, a characteris- 
tic of a FOPT. This phase separation results in a distinct 
abnormal variation of the order parameter. Plausible ar- 
guments indicate that this in turn may lead to a critical 
exponent of 1 or even bigger instead of the normal mean- 
field value 1/2. 

To be specific, we shall speak of a monolayer-like sys- 
tem composed of rodlikc polymer molecules, each of 
which possesses an orientation such that there is a contin- 



uous orientational transition at a certain pressure P. A 
rough feature of the system is that as it is compressed, the 
order parameter 5' of the orientation increases. In addi- 
tion to this orientation transition, however, the molecules 
can change from a gaseous to a liquid state with compres- 
sion. As a result, it is reasonable to assume that only in 
the liquid phase can the orientations of the molecules 
order, since the comparative small average distance of 
the molecules in the liquid renders their interactions 
more effectively. Therefore, our model system should 
be described by the coupling of the first-order gas-liquid 
and the second-order orientation phase transitions. One 
should keep in mind that such a concrete picture is only 
for the sake of convenience: the theory should be gener- 
ally applicable. Thus, the density jump of the gas-liquid 
transition may correspond, for example, to a transition 
from an expanded to a condensed liquid states 0. 

Consider first the pure gas-liquid transition. By taking 
the order parameter C (x) as the area of molecules of one 
mole at the positions around x, the (Gibbs) free-energy 
density for the gas-liquid transition is approximated by 
two pieces of parabolas around the gas and the liquid 
respectively, as 



molar area and Cj* 



/c({C(x)}) = ^«(VC(x))^ 



Ji ' 



for C>Cc 
for C <Cr 



(la) 



with 



Cc being determined by continuity: fg{Cc) — 
fi{Cc), which gives the barrier between the two phases, 
and 

/°({C(x)}) = :| (C(x) - C^f + {P- Po)C{^), 

/°({C(x)}) = ^ (C(x) - C^f + iP- Po)C(x), (lb) 

where Pq the equilibrium transition pressure. We have 
further simplified the free-energy by choosing identical 
positive coefficients /o and k for the two phases. From 
Eqs. Cc (C^ + Cf)/2, independent of P. 

The special forms of the free-energy can be understood 
by their equations of states. For a uniform system at 
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equilibrium, minimizing Eqs. (1) without the gradient 
term, one obtains 

Cg(P)=C°-(P-Po)//o, 
Q{P)^Cf-iP-Po)/fo. (2) 

Consequently, Cg and Ci are the equilibrium molar areas 
at P = Pq, and Cg and Ci those at P. In the extremal 
case at P = 0, the free-energy of the gas fg is lower than 
that of the hquid by PqAC, where AC = C° - C;° is 
the molar area difference between the two phases. Ac- 
cordingly, the gas is stable. As P increases, the differ- 
ence in free-energy between the two phases decreases. At 
P = Pq, the two phases possess equal energies and so can 
coexist. At this pressure, the gas transforms to the liquid 
as the system is further compressed. When P > Pq, the 
liquid is stable and the gas metastable. We presume here 
that nuclei of the liquid phase have already been formed 
possibly by heterogeneous nucleations and so we do not 
consider nucleations anymore below. 

For a nonuniform system at P = Pq, we assume that 
there is only one straight interface dividing the system 
into two parts. Minimizing the total free-energy J fcd^c 
and taking into account proper continuity and boundary 
conditions, one obtains a one-dimensional- like interface 
at a; = separating the gas and the liquid as, 

AC 

Cix < 0, y) = C° + — exp (x/V^) , 
AC 

C(x>0,y)-Cg"-— exp(-a;/V^). (3) 

Note that Eqs. (||) yield correct equilibrium molar ar- 
eas at P = Po for both phases (Eqs. (^) beyond the 
interfacial region of width 2y^. This is obtained at 
the expense of maintaining different internal pressures 
Po - /o\/^AC/2 and Pq + for the hquid and 

the gas, respectively, an artifact stemming from the dis- 
continuous jump in the slopes of the free-energy, Eqs. (1), 
at C = Cc. This jump exerts opposite pressures of the 
magnitude fo^/nAC /2 on both sides of the interface such 
that the total external pressure of the system is again Pq, 
the transition point. However, the artifact will not affect 
our following analyses; we only need a robust linear in- 
terface separating two phases with their respective states 
given by the equilibrium equations, Eqs. (^). 

Next let us turn to the orientation phase transition. 
By defining a local coarse-grained order parameter 5(x) 
as the average orientation per unit area around x, the 
free-energy density for this part alone is assumed to be 

/5({5(x)}) = ^rS^i^) + luSH^) + (V5(x))^ , (4) 

where r = a{PQ — P) and a, u, and k' are positive con- 
stants. We have assumed that S possesses inversion sym- 
metry. Eq. (Q) represents a pressure-induced continu- 
ous phase transition at Pq from a disordered phase with 



S' = to an ordered phase with S = \/ a{P — Po)/u as P 
increases. This gives a critical exponent (3=1/2 defined 
by S" cx (P — Po)'^ 1^. If a system is compressed from an 
area Aq at Pq to ^ at P no matter whether it is in the 
pure gaseous state or in the pure hquid state, Eqs. (||) 
imply that P — Pq cx — A Q , since the total molar 
number of the system is conserved. Therefore the varia- 
tion of S with Aq — A during compression is the standard 
mean-field law with (3 = 1/2. However, these results may 
be radically changed due to the gas-liquid coexistence. 

We now proceed to the coupling of the two phase tran- 
sitions. Within the framework of Ginzburg-Landau the- 
ory, the simplest coupling complying with the symmetry 
requirement of both transitions is given by 

/c5({C(x)}, {5(x)}) = ig (C7(x) - Co) ^^(x), (5) 

where g is the coupling coefficient, and Co a constant. It 
is expected that the smaller the molar area, the greater 
the orientational order. Taking stability into account, 
we have < g < \J2uJq (see Eqs. (0) below). Since we 
have presumed that both transitions take place at the 
same Po, Co may be regarded as an adjustment of the 
transition point, fcs plusing fc and fs constitutes the 
total free energy density /({C}, {S}) of the system. 

The effect of the coupling is most clearly seen from the 
renormalized prefactor of 5'^(x), 

r'-r + 5(C(x)-Co), (6) 

which is position dependent. As C takes different values 
in the gas and the liquid, the orientation phase transition 
can behave differently. Physically, it is clear that in the 
gaseous region, the molecules are far apart. As a result, 
interactions among their orientations are weak and so 
cannot afford order. In the liquid phase, on the other 
hand, they occupy such a small space that they must 
simultaneously order orientationally. Accordingly, in the 
gas, r' > and so S = 0; while in the liquid, r' < 
and so S > 0. This is the most important ansatz whose 
consequences will be explored below. 

The first result of our ansatz is that for a uniform sys- 
tem, the equilibrium solution now becomes 

C;(P) = C,(P)=C,"-(P-Po)//o, 

Sg{P) = 0, 

Q'(P) = C(P)-:^5f(P), 

SfiP) = B{P-Pi)/g, (7) 

where B = g{a + g/fo)/{u — (7^/2/0) is a dimensionless 
constant, and Pi = Pq+ ^ ,' ,f°' is the pressure at which 
a pure liquid would begin to order orientationally. These 
expressions show that at a given pressure, the coupling 
makes the molar area of the liquid reduced because the 
molecules tilt up or straighten on ordering as expected. 
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whereas that of the gas retains their disordered value. 
Meanwhile, there is a finite net orientation in the liquid. 
Using these expressions and Eq. (||), the ansatz that the 
liquid is ordered orientationally and the gas disordered 
can be rephrased as a restriction on P as 



Pi<P<Po 



.9(g° - Co) 
a + g/fo 



P, 



9' 



(8) 



where Pg is the correspondence of Pi for the gas, i.e., 
the pressure at which a pure orientationally disordered 
gas would start to order orientationally had it not trans- 
formed into a liquid. 

Now we determine the pressure where both phases pos- 
sess an equal free-energy, i.e., the transition pressure Pt, 
which must of course be bigger than Pi. To this end, 
it is noted that the orientational ordering of the liquid 
state lowers its free-energy as compared to orientational 
disordering. As a consequence, the relative depth of the 
two free-energy wells is so reduced that now a lower pres- 
sure than Pq can eliminate the difference, since the gas 
remains orientationally disordered. This implies that we 
must have Pi < Pt < P^. This is reasonable because 
otherwise if P; > Pq, the system will turn into the liquid 
state before it becomes ordered. Keeping this in mind, 
one can get from Eqs. (0) and the total free-energy. 



Pt = Pi + 



2(Po - Pi) 



l + ^l+B{Pa-Pl)/{Pg-Pl)' 



(9) 



Note that the transition point of a FOPT is generally 
only the point at which both phases obtain equal free- 
energy. It is not a singularity point of the free-energy. 

Having known Pt , one can determine the constraint on 
Co from the definitions of P/ and Pg and Eq. (|^), i.e., the 
range where Pt is available. First the physical outcome of 
Pi < Pt < Pq implies Co > Cf. Then, on the one hand, if 
Pg ^ P07 or equivalently Co < C^, meaning that a pure 
gas would start to order orientationally beyond Po had it 
not been transformed into liquid, one gets Pi < Pt < Pq 
and so Cf < Co < C°. On the other hand, if Pg < Pq or 
Co > Cg, then Pt must be smaller than Pg. This leads 
to Co < CO + PAC/4, and so C^ < Co < C^ + PAC/4. 
In this case, the coupling pushes Pt to far below Po. 

Here come our main results. As the system is com- 
pressed, the pressure increases. When P reaches Pi, it 
would be possible for the liquid to start to order orienta- 
tionally had it formed. Yet, such liquid embryos possess 
a high energy than the bulk gas and so are doomed to 
dissolve as P is still smaller than Pt. Only at Pt can 
they become stable to further grow (if they are larger 
than the critical nucleus). On the other hand, once the 
liquid nuclei are formed, they must simultaneously be 
orientationally ordered, otherwise they are still not sta- 
ble enough to survive since P is less than Po, at which 
an orientationally disordered (but unstable with respect 
to orientational order) liquid acquires equal energy to 



a disordered gas. Consequently, only those orientation- 
ally ordered liquid regions remains and can live with the 
gas. The system therefore displays a physical picture in 
which an orientationally disordered gas coexists with an 
orientationally ordered liquid with an interface between 
them. The interesting point here is that the continuous 
orientation phase transition exhibits an inhomogeneous 
S'(x) = and S'(x) > phase separation, a characteristic 
of first-order phase transitions. 

As a result of this phase separation, the variation of S 
demonstrates a distinct mode. As the system is further 
compressed, the pressure remains at Pt like an ordinary 
gas-liquid transition. The total area of the system is 
reduced by transforming the gas into the liquid. In so 
doing, the new point here is that the random orienta- 
tions of the molecules get gradually align as they pass 
the interfaces such that the order parameter acquires a 
constant local value Si{Pt) in the formed liquid. The 
physical reason behind is the interactions among the ori- 
entations effect progressively as their average distances 
reduce from the gas to the liquid and then remain at the 
respective constant local values in the two phases. Ac- 
cordingly, as the amount of the liquid increases, the total 
or integrated orientation order of the system increases, 
but Si{Pt) remains constant. It is therefore the expan- 
sion of the total liquid regions with a constant Si{Pt) 
rather than the increase of the local S'(x) itself that con- 
tributes to the increase of the total orientation order of 
the whole system during compression. 

The above analyses are exact and independent of the 
spatial dimensionality D (except Eqs. (|^)). In the follow- 
ing we give plausible arguments to see how the abnormal 
mode of growth of the order parameter is different from 
the normal 1/2-mean-field law. 

At the beginning of the transition, all molecules in the 
system with an area ^0 are in the gaseous state with a 
molar area C'gt (the extra subscript t stands for those 
values at Pt hereafter). If the system is compressed to 
an area A at the constant Pt , the orientationally ordered 
liquid is formed and separated from the orientationally 
disordered gas by interfaces, which is a resultant of the 
gas-liquid and the order-disorder interfaces. We assume 
below that the interfaces have such a thin width (so long 
as K and k' are sufficiently small) that they can be sep- 
arated from the bulk of the liquid and so the mean-field 
approach is valid. Let Ai denote the total area of the 
interfaces, Ci and Si the average molar area and the 
average orientation per unit area in the interfacial re- 
gion, respectively. By noting that the total number of 
molecules is conserved, the total area of the liquid state 
is then 



Ai 



C'gt 



^It ^gt - ^it 



(10) 



a kind of the lever rule with the interface correction. Ac- 
cordingly, the average orientational order parameter of 
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the whole system is given by 

S = {SuAi + S,Ai)/A. (11) 

As pointed out above, within the usual mean-field theory 
of critical phenomena, 5" oc {Aq — A)^ with /3 = 1/2 
independent of the spatial dimensions. We now argue 
how (3 is changed by the distinct growth mode. 

If the system is bisected by a linear interface like 
Eqs. (||), its effective dimensionality is 1. In this case, Ai 
is a constant proportional to the width of the interface. 
Consequently, S is proportional to — A, and so /? = 1. 
For a system of higher dimensions, circular or spherical 
interfaces render At grow with Ai (For D > 2, A means 
(hyper) volume). Assume that the width of the interfaces 
does not change during growth (accordingly, neither do 
Ci and Si). For a hquid region of radius i?, Ai oc R^~^ 
and Ai cx . Accordingly, for not too small size of the 
liquid regions, the bulk contribution Ai is dominant and 
so S is proportional to Aq — A to the leading order, giving 
rise to [3 — 1. For small size, on the other hand, Ai is 
dominant. Accordingly, Ai Ao — A from Eq. ( |l0[ ) and 
so Ai (X {Aq — A)-^/'^-^"^^ to the leading order. However, 
since the factor before Ai in Eq. (^^ is less than one, and 
Sit is larger than Si, also the width of the interfaces may 
change for too small size, the two leading orders with the 
respective exponent D/{D — 1) and 1 of vs {Aq — A) 
may be comparable. Therefore, the critical exponent of 
the order parameter, /3, may appear to be in the range of 
1 and 2 (for D — 2) instead of the normal 1/2. A result of 
this competition at the beginning of the transition (but 
still beyond the critical region) is that the variation of S 
vs Aq — A may become concave rather than the normal 
convex of /9 = 1/2. This result may be relevant to some 
experiments pO) . 

In this paper, we have attempted to provide an over- 
all theoretical understanding of the coupling of phase 
transitions in complex systems by studying a model sys- 
tem consisting of a first-order gas-liquid transition and a 
concomitant second-order orientational phase transition. 
Using a simplified Ginzburg-Landau theory, we are able 
to show analytically that coupling of the two transitions 
may make the continuous orientational phase transition 
exhibit a phase separation of an ordered state from a dis- 
ordered state. In particular, as the system transforms at 
a constant pressure Pt from gas to liquid, the random ori- 
entations of the molecules become ordered continuously 
as they pass through the interfaces between the gas and 
the liquid, and then acquire a constant local value Si{Pt) 
in the formed liquid. Such a phase separation gives rise to 
an abnormal growth law in which the order parameter of 
the orientation transition increases with the expansion 
of the ordered liquid regions during compression. This 
growth law in turn makes the critical exponent P change 
from the 1/2-mean-field value to 1 or even bigger. 

The prerequisite of these results is based on the ansatz 
that only the liquid is orientationally ordered. This re- 



quires that the coupling lowers the transition pressure. 
Critical fluctuations should not change the overall pic- 
ture once Pt is not too close to Pg. Neither should higher 
order couplings since they only renormalize the higher or- 
der terms in the free-energy. 

We emphasize that the theory presented should be ap- 
plicable to other systems though we focus our attention 
only on a monolayer- like system consisting of polymers. 
A salient feature of the complex system under considera- 
tion is that it contains two types of order that can inter- 
play. Another entity besides polymers that may possi- 
bly demonstrate the phenomena discussed here may be, 
for example, electrons, which possess charge and spin. 
When the electrons make up condensed matters as an im- 
portant component, their position order such as charge 
density waves and even the superconducting condensate 
may well compete with their spin order via electromag- 
netic interactions. So one may readily imagine analogous 
pictures to be happened in such systems as well. 

This work was supported by the Postdoctoral Science 
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